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Abstract 

We claim that the non-resonant contribution to non-leptonic charm 
meson decays may not be constant in the phase space of the reaction. 
We argue that this can be relevant for any weak reaction. We discuss 
in detail the decay D + —* K~ir + ir + . 

Non-leptonic charm meson decays have been extensively studied both 
theoretically and experimentally. The high diversity and low multiplicity of 
decay channels provide important information on both weak and strong in- 
teractions. These decays have contributions from resonances in intermediate 
states, as well from the direct non-resonant (NR) decay. The understanding 
of the decay pattern of charm mesons as a whole, and therefore the extrac- 
tion of the decay partial widths for all contributing states, is essential in 
addressing many open problems in charm physics. 

The Dalitz plot analysis |l| is a powerful technique widely used in the 
study of resonance substructures on charmed meson decays. The plot repre- 
sents the phase space of the decay and it is weighted by the squared amplitude 
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of the reaction. Therefore, it contains information on both the kinematics 
and the dynamics. Within this technique, intermediate resonant and non- 
resonant contributions are fitted to get the respective amplitudes and phases. 
The corresponding partial decay widths can then be obtained. 

When experimental data on non-leptonic decays of charm mesons became 
available in the seventies, J. Wiss et al || used the Dalitz plot technique to 
search for the spin of the recently discovered charged D meson. They found 
a result statistically compatible with a flat distribution. Assuming that the 
structure on the Dalitz plot is dominated by the hadronic spin amplitude [§], 
they concluded the D + meson would be a spin particle. 

Subsequently, resonances were found in higher statistics experiments. 
Since then, attention has focused on them and the NR contribution has 
been assumed to be constant. For instance, data on non-leptonic decays of 
the D meson has been fitted |5], ^ 0, |8j using Breit-Wigner functions M to 
represent the various resonances (with the respective angular distribution) 



and a constant function to describe the NR contribution |L0 . 

Although the above parameterization is widely used, a very poor fit has 
been reported [§, |J, suggesting that it may not be adequate to describe 
these decays. These poor results do not improve with higher statistics or 
considering a larger number of resonances ||. Moreover, this problem appears 
in all the D — > Kirn decay channels already measured {D — > K°7r + 7Y~, 
D + -> K°7r + 7T°, D + -> K-n + n+ and D° -> K~n + 7r°) |f| and the worst 
fit is obtained for D + — > K~tv + ii + , where the NR contribution dominates ||. 
(In this case, with 29 degrees of freedom, the \ 2 P er degree of freedom is as 
bad as 3.01.) 

A possible explanation for these discrepancies is the incorrect use of a con- 
stant amplitude for the NR contribution. An incorrect parameterization will 
certainly influence the fit of the resonances and consequently the extracted 
values of amplitudes and phases. As an example, Marklll reported || sig- 
nificant discrepancies on the measurement of the branching ratio (BR) of 
D + — > A > *7r + obtained from the different final states K n°n + and K~ty + ty + . 
Note that while the NR contribution to the first final state is of the order of 
15% of the total partial decay width, in the second it is as large as 80%. 

Here, we claim that NR charm meson decays may contain information 
beyond the simple hadronic amplitude of a spin zero particle decaying into 
three spin zero daughters. Since we are dealing with weak decays, signatures 
of this fundamental interaction can directly appear in the NR amplitude. 
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In weak interactions between quarks and leptons helicity plays an impor- 
tant role. Consequently, one expects a significant dependence of the weak 
amplitudes on the momenta of the interacting particles. Thus, the dynam- 
ics of these reactions vary from point to point of the phase space and the 
significance of this variation depends on the specific physical reaction. 

This should be particularly important in weak decays of charm mesons. 
The large value of the charm quark mass allows for a quasi perturbative 
treatment of QCD. Furthermore, charm quark decays into light quarks and 
this enhances the importance of helicity. For example, we can see the effect 
of weak partonic mechanism responsible for the Cabibbo favored D meson 
decays, i.e. c — > sud, by analysing the decay of r leptons, r — > //P M v T , 
which are essentially similar. This simple example will shed some light on 
the dependence of a weak reaction on its phase space. 

The theoretical Dalitz plot corresponding to the decay r — > /iz/ M z/ r can 
be obtained by taking the well known decay amplitude of pure leptonic 
decays [0]. This decay amplitude can be written as a function of two in- 
variant variables defining a Dalitz plot, e.g., rn 2 p = (p M +Pp m ) 2 and m 2 ^ = 

\M t ^^ Ut \ 2 oc ml UT (m 2 T - m\ VT ) (1) 

where m T is the r mass. 

The dynamics of the reaction has a quadratic dependence on the variable 
m 2 ^. As the Dalitz plot is weighted by | M-t^v^ | 2 , equation ([I]) shows that 
a Dalitz plot of a pure weak decay has indeed significant variations along the 
phase space. 

Obviously, due to the hadronization process of the partons after their 
weak interaction, the result of the previous example cannot be simply trans- 
lated into hadronic decays. In the latter case, one has to take into account 
non-perturbative QCD effects involved in the final hadronic state formation. 
In order to make an estimate of the effect of the dynamics in the Dalitz plot, 
we use an approximate method to describe hadronic decays. The method is 
based on both the factorization technique []13] and an effective Hamiltonian 



, [T5|] for the partonic interaction and has been successfully used to describe 



heavy meson decays 13, 16 



As we are interested in the NR contributions, we analyse the channel 
D + — > K~n + 7Y + , which has a very large NR branching ratio, as mentioned 
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above. The effective Hamiltonian for the weak vertex c — > sud is [H] , |T3|: 



Gp 

Tl e ff = ( — j=) cos 2 ^cfo-i : (sc)(ud) : +a 2 : (sti)(Sc) :] (2) 

where (gg 7 ) is a short-hand notation for g7 M (l~ 7s)?'- The coefficients a% and 
02 characterize the contribution of the effective charged and neutral currents 
respectively, which include short-distance QCD effects. Their values have 
been fitted in the case of charm meson two-body decays (see for example 
reference [fUjl ). The diagrams contributing to the decay D + — > K~n + n + are 
shown in Figure (). Using factorization we obtain the following decomposi- 
tion for the hadronic amplitude 

M D +^ K - 7T + 7T + = (—/=) cos 2 9 c [ai(K~7cf |sc|-D + )(7r^|"u<i|0) 

+a 2 (K-TT+\sd\0}(7r+\uc\D + ) + (7^ <-> tt+)] . (3) 

Let us first discuss the term driven by ai, i.e, the one of Figure (.a). The 
most general form to decompose the first matrix element can be written in 
terms of four form factors 0]. Using the parameterization of reference [fTS], 
we can write: 

(K-ir+\sc\D+) = A^F 1 + A%F 2 + iV^F z + A%F 4 , (4) 

where 

,Q ■ (pk + Pd) 



The terms proportional to F±, F 2 and F4 originate from the axial vector 
part of the matrix element whereas the one proportional to F3 originates 
from the vector part; the terms proportional to F\, F 2 and F$ correspond 
to spin 1 and F4 to spin 0. The four form factors depend on three variables 
m i = (Pfc +Ptti) 2 ? m \ — (Pk + Pn 2 ) 2 an d Q 2 which is a constant (m 2 ) in this 
case. 
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The second matrix element in equation (|3]) has the well known form 

(n+\ud\0) = iUp% 2 . (5) 

The only contributing term in equation (||) after multiplying it by equa- 
tion flU), is the axial spin term, i.e., 

(K-7rf\sc\D + )(7r+\ud\0) = (k 2 „ F A ) {if vPv ») = %j^m\F± . (6) 

To find the contribution of Figure (.b), one can use the well known 
expressions jl9[ 



(ir^\uc\D + ) 



{Pd + V 



+ 



m D — m; 



and 

(K-(p K )7rf\sd\0) 



(7if\sd\K + (-p K )) 



,-pK+Pm, 



m 



K 



m 



~(-PK -P* 



U(q 2 



+^ L V^(-PK-p^rfo(q 2 ) 



In the equations above, q 2 = (po — Pn 2 ) 2 — (~Pk — PttJ 2 while the 
functions F^ n (q 2 ) (corresponding to a current of spin parity J p ), f+(q 2 ) and 
fo{q 2 ) are form factors. We return to them later. 

We then find for the second contribution in equation ([$]), 

(7r + \uc\D + )(K-n + \sd\0) = F^{m\) j ' + {m\) (m 2 D + m 2 K + 2m\ - 2m 2 2 - m\) 

I [pl-Z^w /_2\ p0+/„2w / 2\l (j^JT m2 )( m K ~ m D 
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mf 



+ (ml <-> ml) 



(7) 



where we have explicitly introduced the Dalitz plot variables m\ and m\ 
defined above. 
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The contribution of diagram (.a), given by equation (|6]) is proportional 
to f^m 2 . Thus, unless the form factor F 4 is unacceptably large (F 4 ~ 10 3 ), 
we can safely neglect this contribution in favor of that of diagram (.b), given 
by equation (0) which contains m 2 D . As an aside, it is possible that the NR 
part of the decay D + — > K~it + it + is large precisely because the contribution 
of diagram (.b) is not small. 

The NR contribution to the amplitude of the decay D + — > K^n + n + can 
thus be simply written replacing equation (|7|) in (|3|), neglecting the contri- 
bution of Figure (.a). The final expression thus depends on the effective 
coefficient 02 and the four form factors. The two Dit form factors F^lq 2 ), 



have well established expressions ||15|| : 

where M Dnjl - = 2.01 GeV and M B7r + = 2.2 GeV. They have been suc- 
cessfully used in the kinematic range we are considering here. The poles 
lie outside our kinematic region. The Kn form factors, f+{q 2 ) and fo{q 2 ), 
can be extracted from the semi-leptonic decays K — ► nlu, with I = e, \x. 
Nevertheless, it is not clear that the usual parameterization ]T2" | 



f + {q 2 ) = / + (0) ( 1 + A + 41 , fo(q 2 ) = /o(0) I 1 ■+ A (1 -^ ) (9) 



m 2 




is valid in the whole kinematic region of our reaction. In equation @, /+(0) = 
/o(0) = 1 and the other coefficients have been measured to bepOf: A + ~ 0.03 



independent of the measured channel, whereas the value of Ao depends on 
the decay: A ~ for K~ — > Tx^yTv and A ~ 0.025 for K° — > 7r + fi~v. 

In order to check the validity of this calculation scheme, we have evaluated 
the NR partial decay width T(D + — > K^tv + tv + ) nr using the expressions 
above. With A = and the value of a 2 extracted from two body decay fl4 



we find a branching ratio (BR) of 9% which is close to the reported exper- 



imental value^OJ 7.3 ± 1.4% obtained by fitting the NR contribution to a 
constant. We studied the stability of this result under the change of the 
parameters A + and Ao : if we take the various values extracted from different 
channels we find that the BR varies less than 30%. Even assuming constant 
form factors (A+ = Aq = 0), the BR remains of the same order of magnitude. 
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Figure () shows the Dalitz plot for the NR contribution to the decay 
D + — > K~n + n + as a function of the variables ml and m\. It has been gen- 
erated by Monte Carlo simulation with a weight proportional to the square 
of the amplitude in equation @, using equation (|7|). We have considered the 
same central value of the parameters as above. As one can see from equation 
fl7|) and Figure (), according to this calculation the matrix element describing 
the dynamics of the NR contribution to the decay D + — ► K~-k + 'k + signifi- 
cantly varies along the phase space of the reaction. Its shape remains almost 
the same for other values of the parameters of the Kn form factor. This is 
still valid even if we take the four form factors as constants. 

However, the result presented in Figure () has been obtained using an ap- 
proximate method. Non-perturbative effects, present in this decay through 
the exchange of soft gluons, or final state interactions could change the struc- 
ture shown in this figure. In the extreme case where non-perturbative effects 
completely dominate the decay, the structure will be washed out because of 
the dispersive nature of these effects, therefore obtaining the flat contribu- 
tion predicted by the pure hadronic decay of a zero spin particle. Comparison 
between the distribution shown in Figure () and experimental data will thus 
be a test for the validity of the factorization method. 

In summary, we have shown that the natural parameterization for the 
non-resonant part of charm decays - based in the spin amplitude of the 
hadronic decay - could significantly change due to the fundamental weak 
interaction between quarks. The appearance of these structures in the plot 
could be responsible for the problems of the various experimental teams with 
the convergence of their fits. To clarify this point, it is important in future 
analyses to use a parameterization for the non-resonant contribution going 
beyond the simple constant. 
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Figure 1: The two diagrams contributing to the decay D + 
according to the effective Hamiltonian of equation (2). 
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Figure 2: The Dalitz plot of the decay D + — > i^7r + 7r + , weighted by 
|AlD+-*ft:-7r+7r+| 2 as in equations (3) and (7), generated via Monte Carlo. 
The Dalitz plot variables are m\ = (p K + p^) 2 and ml = (p^ +Pvr 2 ) 2 - 
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